clearvars -except Globaloption option 
label = '-WithDebt';

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% This code reads in the VAR dynamics and produce the Campbell-Shiller decomposition in the paper
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Load data

loaddata_4CIjoint_annual_debt;
ttime=[startdate:enddate];

output_rp = 0.03;

%% Estimation
run ../../tools/cs_estimation;

%% steady-state upper bound calculation
upper = exp(pxbar)*(mean(taxrevgdp-spendgdp));
save MAT/debtinVAR_UK_pre1945.mat s gdebt upper

